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ABSTRACT 

We present infrared photometric observations of the Be binary system 5 Scor- 
pii obtained in 2006. The J,H and K magnitudes are the same within the errors 
compared to observations taken 10 months earher. We derive the infrared ex- 
cess from the observation and compare this to the color excess predicted by a 
radiative equilibrium model of the primary star and its circumstellar disk. We 
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use a non-LTE computational code to model the gaseous envelope concentrated 
in the star's equatorial plane and calculate the expected spectral energy distri- 
bution and Ha emission profile of the star with its circumstellar disk. Using 
the observed infrared excess of 5 Sco, as well as Ha spectroscopy bracketing the 
IR observations in time, we place constraints on the radial density distribution 
in the circumstellar disk. Because the disk exhibits variability in its density 
distribution, this work will be helpful in understanding its dynamics. 

Subject headings: Stars, Astronomical Techniques 



1. Introduction 



Be stars are massive, rapidly rotating stars that exhibit distinct characteristics due to 
the presence of a thin disk of gas concentrated in the stars' equatorial plane. The classical 
definition of a Be star is that of a main or near-main sequence B star whose s pectrum has, 
or had at some point in time, conspicuous Balmer emission lines (j Collins! 1 198 71 ). In addition 
to this prominent emission line spectrum, important observational features of these rapid 
rotators include a partial polarization of radiated light and a continuum excess in the radio 
and infrared (IR) spectrum. These distinct observational characteristics of many Be stars 
are variable on timescales of weeks to decades, due primarily to changes in the physical 
conditions in the circumstellar disk. The preceding description limits classical Be stars to 
intermediate mass stars surrounded by gaseous envelopes and excludes related stars such as 
the Herbig Ae/Be stars which are pre-main sequence objects surrounded by dusty disks. 

The contribution of radiation from the circumstellar gas to the Be stellar continuum 
is significant. Observations of Be stars at IR wavelengths show a much brighter continuum 
spectra than would be expected from typical stellar photospheric emission. The explanation 
for this IR excess is that the continuum of the cent ral star in the IR is dominated by free- free 
and free-bound emission from the equatorial disk (IGehrz et al.lll974j ). 



While evidence suggests that the gaseous, equatorial disk is formed from material re- 
leased from the star, the physical proces ses that govern the formation of the disk are not 
well understood (jPorter k. Riviniusll2003l ). and are the subject of persistent study. The first 
viable models predi cted an equa torial mass-loss disk triggered by the near-critical rotation 
of the central star (jStruvelll93ll ). While there is general agreement that the rapid rotation 
of the central star is fundamental to the disk formation mechan ism, its proximity to the 
critical value remains a contentious topic (ITownsend et al.l 120041 ). Other models of stellar 
mass ejection suggest that the gas outflow could be due to radiatively driven winds, possibly 
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triggered by magneti c fields (jPorterl 119971 ) . or due to pulsations in the stellar ph otosphere 



Rivinius et al.ll200ll ). Recently, the viscous decretion disk model pioneered by iLee et al. 



( Il99ll ) has successfully reproduced several key observables, including the IR excess, from Be 



star disks (jPorterl Il999l ) 



The Be star 6 Sco is of particular interest to those studying the mechanisms of disk 
formation. In the surnmer o f 2000, the star began exhibiting line emission in its radiated 
spectrum (jOtero et al.l 120011 ) . Spectroscopic monitoring of the star over the following two 
years confirmed a continuous strengthening of the emissi on lines due to the steady develop- 
ment of a circumstellar disk (IMiroshnichenko et al.ll2003l ). In addition, the bri ghtness of the 



star has varied considerably, increasing by as much as 0.7 mag in the V band (iGandet et al. 



2OO2I ) . Because 6 Sco presents a unique opportunity of observing a Be star in an active disk 
formation phase, it is of particular interest for detailed astrophysical studies. 

In this paper, we present unique infrared observations of the star 6 Sco and we compare 
the observed IR excess to that predicted by our theoretical models. Our theoretical models 
are generated by a non-LTE code that enforces radiative equilibrium by a self-consistent 
disk thermal structure and includes a solar-type chem ical coraposit i on. For this investi- 
gation, following p revious work to match the IR excess (IWaterslll986l : ICote fc Waters! 119871 : 
Waters et al.lll987l ). we have adopted a simple density distribution which falls off as a power- 
law with increasing distance from the central star. Our models are also constrained by com- 
pari son of predicted Ha profiles with contemporaneous observations. It should be noted 
that ICarciofi et al.l (120061 ) reported that 6 Sco's density distribution during this active stage 
of disk growth seems to deviate from a simple power law. However, as pointed out by 
Okazakil (120071 ) . in order to assess the validity of the viscous decretion model it is important 
to carefully monitor the IR excess of Be stars to determine if the density index changes with 
time. Ascertaining a reasonable average value of the index at the time of IR observations, 
in conjunction with other work, will be an important test of the viscous disk model. In the 
following section, we discuss our photometric observations of 6 Sco. In section 3, we present 
our models for the circumstellar disk around the studied star. Lastly, we summarize the 
results of our work in section 4 and discuss the implications of our findings. 



Observations 



S Sco is a relatively bright (K=2.61 ± 0.07 mag.. lNeugebauer &: Leightonlll969l ) star and 
presents difficulties for observations with anything but small aperture telescopes with modern 
array detectors. The telescope and detector used were the 61" (1.55 m) Kaj Strand astro- 
metric reflector and the Astrocam infrared array camera. The Strand reflector is an F/9.8 
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folde d Newtonian with a plate scale of 13.5"mm ^ (jStrandlll97ll ). Astrocam (IFischer et al. 



20031 ) contains a 1024 x 1024 InSb array with a minimum readout time with double corre- 
lated sampling for the full array of 0.19 seconds. Astrocam pixels are 27 n square providing 
0.37 arcsec sampling with 1:1 reimaging of the focal plane. Observations were made on 2006 
May 10, 11 and 12. The first night was used as an engineering test run to determine the 
best strategy for observing 6 Sco. Of the following nights only 2006 May 12 was deemed 
photometric. Initial tests showed that even with the minimum readout time Astrocam im- 
ages of 6 Sco were saturated. A +10 mm defocus was applied which produced the necessary 
donut shaped images (Fig. [1]) to avoid saturation with 0.20 sec exposures and 50 co-adds in 
all filters (J,H,K). The full width of the donuts was measured in several images and found 
to be 38 pixels, or 1.026 mm. This is in excellent agreement with the geometrical optics 
expectation of 1.0204 mm for the Strand telescope. The central spot visible in Fig. [T] is 
known as Poisson's spot and can be understood from diffraction (see Hecht and Zajac 1979 
for the source of this ironic name). All data and calibration images were obtained with the 
same +10mm defocus. We elected to observe both S Sco and the photometric s tandards with 



this sa me defocus to simplify our analysis. Standard stars were chosen from iGuetter et al. 



(120031 ) and are listed in Tabled] with their K magnitudes, J — H and H — K colors, exposure 



times and number of co-adds which were the same for all filters. 



Data reduction and image processing were performed with IRAF (|Todvlll993). Ind ivid- 



ual frames were linearized following a procedure described by iLuginbuhl et al.l (Il998l ) and 
then flat-fielded. All stars were observed with a 4 position dither pattern and these images 
were used to construct sky images for subtraction. A bad pixel mask was created and used 
to correct the images. After this process was performed each star was identified in the image 
and coordinate lists prepared. Aperture photometry was performed to produce curves of 
growth to determine the best parameters for aperture photometry. At this point we noticed 
that the 4 dithered images generally produced offset pairs of curves of growth. This offset 
was at the ~ 5% level and we attribute this to the different amplifiers used on Astrocam. 
This offset was always visible between the top and bottom halves of the array. We increased 
the counts in pixels below the camera midline by 5% to compensate for this offset; Fig. [2] 
shows the curves of growth before and after the 5% correction was applied. For an aperture 
radius greater than 50 pixels we generally saw a divergence in the curves of growth and we 
therefore selected 50 pixels as the photometric aperture. We attribute this divergence to 
the difficulty in performing an accurate sky subtraction across large areas of the array and 
believe this ultimately limits our photometry. 

A final factor to contend with were some cosmetically poor areas of the array. These 
areas had very different noise characteristics from the rest of the array. There was a significant 
positive bias in the noise which artificially elevated the count levels. The exact amount of 
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the increase varied as a function of position in these roughly rectangular shaped areas. As 
this was very difficult to calibrate out we elected to disregard images for which any portion of 
the photometric aperture and sky annulus fell within these areas. Stars could be artificially 
brightened by up to 0.2 magnitudes if they fell within these areas. This restriction eliminated 
14 out of 96 standard star frames and 2 out of 48 6 Sco frames. 

The aperture photometry results were exported to custom software written in IDL. Fits 
for the extinction were made in each filter and are shown in Fig. [21 The initial formal 
errors were determined by adding in quadrature t he error in the instru mental magnitudes 
determined from photon statistics and those from iGuetter et al.l (120031 ) . These initial fits 
yielded large reduced x^- We elected to increase the error bars by a factor of the square root 
of the reduced and these errors are shown in Fig. [31 We believe this increase in the errors is 
needed to properly account for the errors in estimating the sky subtraction across large areas 
of the array at the four dither positions. We corrected the standard star results for extinction 
and then tested for the presence of any significant color terms. Fig, [j] shows t he de-extincted 
magnitude differences plotted with the J—H colors from lGuetter et al.l (120031 ). No significant 
color terms are apparent by eye and the color term slopes with errors are consistent with 
zero (Table [2]). The reddest standards could not be observed in K after a computer disk 
unexpectedly reached capacity late in the night of 2006 May 12. We do not feel that the 
slight positive correlation for K evident in Fig. [H justifies the inclusion of a K color term. We 
expected a negligible color term as the filters were originally selected by considering their 
transmission and the quantum efficiency of the Astrocam detector to minimize any color 
term (A. Henden, private communication). Table [2] summarizes the extinction and color 
term results, and the increase to the photometric errors applie d in Fig. [3l Our ext inction 
coefficients are in reasonable agreement to those determined by iGuetter et al.l (120031 ). 



The J data in Fig. [3] exhibit a larger scatter than the other filters. This could be due to 
a variable water vapor content. Although water is not the only source of o pacity it is a sig- 
nificant effect and is more pronou nced in the J filter compared to H and K (iManduca &: Bell 
19791 : IVolk. Clark fc Milondll989l ). Fortunately a ground based Global Positioning System 
(GPS) meteorology station is located on the Camp Navajo Army Depot 8.3 km WNW from 
the USNO-Flagstaff Station . This station (station code FSTl) is part of the National 
Oceanic and Atmospheric Administration GPS Integrated Precipi table Water (IPW) Vapo r 
(NOAA GPS-IPW) network. The system is described in detail in lWolfe k Gutmanl J200o[ l. 
The IPW content in centimeters over the site is recorded every 30 minutes, and the data are 
available via the world wide wet0. 



^Observations courtesy of the NOAA Earth Systems Research Laboratory: 



http://gpsmet.noaa.gov 
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Our IR observations of 5 Sco and the standards were undertaken between 0557 UT and 
0853 UT. During this time the IPW content decreased from 7 to 6 mm (Fig. [5]) at the Camp 
Navajo station. While it is beyond the scope of this work to quantify the results of this 
decrease we can qualitatively expect a smooth decrease to the opacity over time based on 
the GPS IPW monitoring. To explore the effects of a variable water vapor content we show 
the de-extincted J magnitudes as a function of UT time in Fig. |6l We also display in this 
figure a solid line which shows the average of the de-extincted magnitudes as a guide to 
identify any correlation with the water vapor content. It is also useful to consider Fig. E]in 
conjunction with the airmass of the observations which we show in Fig. [71 For example if all 
the high airmass standard star observations occurred early or late in the night they may be 
biased by a variable opacity and significantly effect the linear extinction fits. 

Considering Fig. [H] a slight brightening may occur from 6 to 8 hours UT, qualitatively 
consistent with a decreasing water vapor content. Furthermore, our 2 observations near an 
airmass of 2 are closely spaced in time but after correcting for extinction they differ in Fig. [6] 
by approximately 0.05 magnitudes. The scatter in the standard star magnitudes in the 
individual sets of dithered images would, if due to a variable opacity, suggest a rather high 
frequency variable water vapor signal. Although the time sampling in Fig. [S] might preclude 
detection of such a signal some large 1 mm IPW swings should have been apparent even 
with the 30 minute sampling interval. While the general trend in Fig. [6] may be indicative 
of a changing water vapor content we cannot conclusively explain the scatter within the 
individual sets of dithered images except to invoke a problem in the sky subtraction. This 
we noted earlier when we saw the divergence in the curves of growth beyond a radius of 
50 pixels (Fig. [2]). Knowledge of the water vapor content ca n be used wit h a model of 



total atmospheric extinction to calibrate out extinction effects ( lAngiondll989l ) with the aid 
of a radiometer. The proliferation of GPS-IPW monitoring sites and their provision near 
astronomical o bservatories could be used to calibrate out effects of atmospheric extinction 



( lGutmanll2007l ) and it will be interesting to explore this possibility in the future. 



Using only the extinction terms in Table [2] we determined the magnitudes of 6 Sco for 
46 images. Fig. [8] shows these individual measurements which within the errors appear flat 
as a function of airmass. The errors are dominated by the errors in the extinction terms 
due to the large counts in the S Sco frames. We averaged the data and reduced the errors 
by and found formal errors of 2% or less in t h e indi vidual magnitudes. Examination of 



Columns 9 and 10 of Table 1 from iGuetter et al.l (120031 ) allows us to establish that a lower 



level floor of 3% is a more probable uncertainty for the individual magnitudes based on their 
many nights of observations at the Flagstaff Station. We therefore adopt this as our final 
uncertainty in the individual magnitudes, and propagate this in quadrature to produce the 
uncertainties in the colors of 6 Sco. Table [3] lists the magnitudes and colors with this final 
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uncertainty es t imate . The magnitudes in Table 3 are identical within the errors to those of 
Carciofi et al.l (120061 ) observed 10 months earlier. 



Theoretical Be Star Models 



Several models of the thermal structure of Be star disks exist that can reproduce the 
spectral energy distribution of a star with its surrounding material. The theoretical models 
produced here were constructed with BEDISK, a non-LTE code that predicts the thermal 
structure of the disk using a re alistic chemical comp osition. The principal features of BE- 
DISK are discussed in detail in ISigut fc Joned (120071 ). We present here a brief overview of 
the code relevant to this investigation. 



3.1. Modeling a Be Star with BEDISK 

The stellar radiation that reaches any particular point in the disk of the Be star depends 
on the emissivity and the opacity of the gas present throughout the entire disk. Thus, 
the photoionizing radiation field at a particular location in the disk is dependent on the 
density, temperature and atomic level populations at other locations in the disk. For a 
given photoionizing radiation field at each location in the circumstellar disk, the code solves 
the equations of statistical equilibrium (non local thermodynamic equilibrium) for the level 
populations for each atom and each ion included in the computation. Having produced the 
atomic level populations, the rates of heating and cooling due to the various atomic processes 
in the material can be computed at each disk location. Since the disk is assumed to be in 
radiative equilibrium, the temperature at each disk location can be iteratively obtained by 
balancing the rates of energy gain and energy loss. As the temperature profile is adjusted, 
the density profile is also adjusted in order to maintain vertical, hydrostatic equilibrium 
throughout the disk. 

It is assumed that the radial density distribution of the disk falls off as a power law in 
the equatorial plane described as: 

p{R) = po{R/R.)-'' (1) 

where po represents the density of the circumstellar material closest to the central star and 
R > R^. This p o wer law parameterizati on of the circumste l lar di sk follows the prescription 
of lWatersI Jl986h . lCote fc WatersI Jl987h . and lWaters elall Jl987h . who constructed models 
to match the IR excess of Be stars. A second assumption is that the circumstellar gas is in 
vertical isothermal equilibrium perpendicular to the plane of the disk. The code accepts a 
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user defined set of atomic models wliicli contain a list of the energy levels and transitions for 
each atom and ion included in the disk model. The BEDISK code incorporates 425 levels 
and over 1 000 radiative and co Uisional transitions for nine elements and thirty-one ions (see 
Table 1 in lsigut fc Jo^boOTh . 



3.2. Ha Emission Line 



The JHK photometric observations of 5 Sco were obtained on 2006 May 12. To measure 
the shape and overall strength of the Ha emission line close in time to the IR observations, 
we have obtained spectra of 5 Sco on 2006 April 19 and 2006 May 22. The spectroscopic 
observations were obtained using a fiber-fed echelle spectrograph connected to the John 
S. Hall 1.1 m telescope at the Lowell Observatory. The spectra in the Ha line region were 
obtained at a resolving power of ~10,000 and with pixel spacing of 0.022 nm (resulting in 
3 pixels per resolution element). The signal-to-noise ratio (SNR) in the continuum level 
near Ha line is approximately 100 for both spectra, and therefore the SNR of the emission 
line is around a few hundred. The spectroscopic ob servations were reduced using standard 
reduction routines written in (IHall et al.lll994j ). and the resulting Ha profiles of 5 Sco 



obtained on two different nights in 2006 are shown in Figure 9. We have not corrected 
for the effects of telluric lines, which can be seen in the Figure as shallow absorption lines, 
however, we have reduced their effect on the equivalent width (EW) measurement of the Ha 
emission line by excluding the contributions from regions that go below the continuum level 
near the emission line. The resulting EW measurements for 5 Sco are -20.0 ± 0.06 A on 
2006 April 19 and -13.8 ± 0.04 A on 2006 May 22. The uncertainties associated with our 
EW measures are estimated based on the limitation of the continuum level determination 
(which we estimate to b e arou nd 3%) and are not associated with the SNR of the continuum. 
VoUmann fc Eversberd (120061 ) discuss the error analysis in equivalent widths in the case of 
photon noise limited case, and our uncertainties are equivalent to their analysis when the 
uncertainty associated with the continuum estimation dominates (i.e., when the second term 
in their eq. 5 is the main contributor to the uncertainty). 

The model Ha profiles were determined by formally solving the transfer equation along 
a set of rays through the system viewed at an angle of i relative to the line-of-sight {i = 0° 
is pole-on, and i = 90° is equator-on). We adopt a value of i = 30° to c o mput e the 
profiles. This i is in agreem ent with the values quoted by both ICarciofi et al.l (120061 ) and 
Miroshnichenko et al.l (l200ll ). The disk is assumed to be in pure Keplerian rotation, and the 



^Interactive Data Language of RSI, ITT Industries, Boulder, CO 
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equatorial rotational velocity of the star was chosen so that the measured v sin i of 150 km/s 
( Miroshnichenko et al.|[200ll ) of the star was consistent with i = 30°. For rays terminating 
on the stellar surface, a photosp heric Ha line profile is c ompu ted using the synthe code and 
the Stark Broadening routines of lBarklem &: Piskunovl (120031 ): these profiles are based on the 
hydrogen pop ulations c ompu ted for the appropriate LTE, line-blanketed model atmosphere 
adopted from iKuruczl (119931 ) . In the disk, the hydrogen emissivity and opacity followed 
directly from the hydrogen level populations determined using the self-consistent thermal 
structure. As a result, the full dependence of the Ha hne emissivity and opacity on the 
physical conditions through out the disk were retained. The Stark broadening routines of 



Barklem &: Piskunovl (120031 ) were also used to compute the local Ha line profile throughout 



the disk. This is important as these routines can handle high disk densities for which the 
assumption of a simple Gaussian or Voigt profile for the Ha line would be invalid. The Ha 
line profile was obtained by summing over the rays weighted by their projected area on the 
sky. Finally, the profile was convolved with a Gaussian having FWHM of 0.656 A so that 
the predicted profile matched the resolving power of 10"^ of the spectroscopic observations. 



3.3. Modeling 5 Scorpii 



We adopted the stellar parameters chosen by ICarciofi et al.l (120061 ) in their study of the 
properties of 6 Sco's circumstellar disk. The radius and mass of the star are 7.0 Rq and 14.0 
Mq, respectively, and the star's effective temperature is 27 000 K. We constructed models for 
seven values of n (n = 1.5, 2.0, 2.5, 3.0, 3.5, 4.0 and 4.5) with varying values of po beginning 
with 1.0 X 10^^^ g cm~^ and increasing the density up to 5.0 x 10^^° g cm~^. We then 
integrated the predicted spectral energy distribution of the stellar models over the J, H and 
K bandwidths in order to calculate the theoretical J-H and H-K colors. The results of our 
work are plotted in the color-color diagram illustrated in Fig. [TOl 

In order to determine our best-fit model, we computed theoretical Ha profiles for a 
subset of models with synthetic colours to the observed value. We compared the equivalent 
widths and line shapes of these predictions to Ha spectra bracketing the IR observations 
in time. These models corresponded to a range of n from 3.5 to 4.5 and po from 5.0 x 
10~^^ to 5.0 X 10~^° g cm~'^. The models with the larger values of n have a greater drop 
in density with radial distance from the star and therefore correspond to the larger values 
of Po quoted. Within this subset, the theoretical equivalent widths range from a minimum 
of 7.1 A to a maximum of 18.8 A for i = 30°. Changing the value of i by ± ~ 10° for 
this subset of models changes the equivalent width by less than ~ 1 A. However, changes 
in i do significantly alter the shape of the profile. Furthermore, we produced our emission 
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line profiles with the inclusion of a hydrostatic subroutine that enables our code to produce 
fully self-consistent disk densities and temperatures. Fig. [9] compares the best-fit predicted 
Ha line profile to the observed lines. Note that the observed line changes dramatically in 
strength and shape between the two bracketing observations taken about one month apart. 
This indicates that there are ongoing changes in the disk density structure during this short 
time period. 

As can be seen from the figure, the computed line profile for the i = 30° model is too 
wide in the wings as compared to the observations. To attempt to reduce the width of 
the line, we have computed the Ha line profile for a slightly smaller inclination of the disk 
relative to the line of sight {i = 20°) and also added an evacuated annulus between 1.5 and 
4.0 stellar radii which had a density of 0.01 compared to our single power-law in this region. 
We acknowledge th at the smaller inclinat i on of i = 20° seems inconsistent with the orbital 
solution for 6 Sco (IMiroshnichenko et al.l 1200 ll ) which gives the orientation of the orbital 
plane relative to the line of sight as 38° ± 5°, assuming that the disk lies in the orbital plane 
of the system. The smaller inclination angle and the addition of the evacuated annulus do 
reduce the width of the line with the i = 20° hole model coming closest to the observations. 
Nevertheless, we are not proposing the more complex hole model as the best fit as there 
are a number of parameters in this model (location of annulus and density contrast in the 
annulus) which would need additional observational constraints. In addition, we note that 
the observed lines show a definite asymmetry which cannot be modeled by our axisymmetric 
calculations. Further study of the line asymmetry could reveal a dynamic eiiyelope density 
structure such as a one-armed oscillation in the equatorial disk (jVakili et al.lll998l ). In the 
future, we plan on incorporating density perturbation into our disk models so that we can 
replicate the asymmetric line profiles. 



In comparing our work with the study on S Sco presented by ICarciofi et al.l (120061 ). we 
emphasize our agreement with some of their results. First, we note that the disk temperature 
profile for our best-fit parameters (Fig. [TTj) shows reasonable agreement with the average 
temperature distribution of their best-fitting model. Both theoretical models predict a cool 
dense region along the equatorial plane near the star. We obtain a minimum temperat ure 
of 6700 K in agreement with the minimum temperature quoted by ICarciofi et al.l (120061 ) of 
7000 K. The average overall tempe rature we obtain is 12000 K, slightly less than the value 
of 16000K of ICarciofi et al.l (120061 ). We plan on performing a detailed comparison of our 
computational code with Carciofi et al.'s non-LTE Monte Carlo radiation transfer code in 
a future work. Furthermore, our predicted colours of the star from our set o f models cover 



the ra nge of values in which we find the photometric observations presented in ICarciofi et al. 
( 120061 ) (Fig. [T2l) . The observations, taken over a five year period during which the disk 



surrounding 6 Sco was developing, suggest that our radial drop-off model for the structure 
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of the disk is a reasonable approximation for the inner portion of the disk where Ha and the 
near IR continuum form. Notice that our theoretical models constructed with n = 4.0 for a 
range of po fall within the locus of the Carciofi et al. (2006) photometric observations. 



3.4. Binary Component 



In order to ensure that our predicted IR excess for the disk of 5 Sco is an accurate 
measure, we took into consideration the binarity of the system by assessing the impact of 
the secondary on the observed infrared continuum. We set the sp ectral type of th e secondary 
star through its optical magnitude difference from the primary. iBeddingl ( 19931 ) found that 
the secondary was around 1.5 mag fainter than the diskless BO. 2 IV primary. Our estimate 
of a B3 spectral type for the secondary is supported by comp aring the estim ated mass of 
the obje ct (the mass ratio of the system is 1.5 - 2 according to IBeddingl Il993l ) to the values 
listed in lAUenl (119991 ). Assuming this spectral type for the secondary, we re-calculated our 
predicted colours by including the SED from the secondary star. Fig. [12] shows the spectral 
energy distribution of the secondary, the primary and the primary with its circumstellar 
disk. 

The inclusion of the secondary star did not change the general range of the power-density 
exponents. For each value of ra, we found that this approximation did not shift the value of 
our best-fit disk-base density by more than 0.03 mag in H-K and 0.02 mag in J-H. While the 
binarity of the system does affect the observed infrared continuum, our investigation shows 
that the influence of the secondary is not considerable enough to significantly alter the range 
of parameters that match the infrared excess observation. 



4. Conclusion 



In summary, we found that 5 Sco's J, H and K mag r iitude s are identical within the errors 
to those observed 10 months earlier by ICarciofi et al.l (120061 ). We modeled 5 Sco's circum- 
stellar disk and compared the theoretical infrared excess to photometric observations of the 
excess. These observations, combined with contemporaenous spect roscopic obseryatiori s, are 
particularly important for studying Be stars in their active phases ( ICarciofi et al.ll2006l ). We 
used a non-LTE computational code to model the thin equatorial disk of gas surrounding 5 
Sco, employing a simple power-law for the density distribution of disk material. Although 
the stellar activity suggests a more complicated density distribution in the circumstellar disk, 
monitoring the infrared excess and determining the average value of the power-law index re- 
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mains an important part of studying 6 Sco (jOkazakill2007l ). With this in mind, we calculated 
the spectral energy distribution of the star and its disk and we derived the theoretical infrared 
color excesses. We compared these predictions to our infrared photometric observations of 
6 Sco and further constrained the disk parameters by modeling the Ha emission line and 
comparing it to spectra taken before and after the time of the photometric observation. Our 
best fit parameters for the circumstellar disk were found to be po = 5.0 x 10~^° g cm^^ with 
a power-law density exponent of n = 4.0. Overall, we found acceptable parameters within 
the range for n from 3.5 to 4.5 and po from 5.0 x 10^^^ to 5.0 x lO"^'' g cm~^. Since 6 
Sco's disk changes dynamically, our analysis places constraints on the disk at the time of the 
observations. It will be interesting to see if this value for the power-law index changes with 
time which will provide evidence to assess the applicability of the viscous disk model. Re- 
cently several investigations h ave suggested speci fic detailed features in density distributions 
in Be star disks. For example, IZorec et al.l (120071 ) find that the power-law index is less than 
one within the first 3 stell ar radii. This i mplies that the density distribution near the star 
is constant or increasing. iRivinius et al.l (120011 ) proposed that Be star disks are comp osed 
of a series of r i ngs. T his concept has been studied recently by by iMeilland et al.l (120061 ) and 



Carciofi et al.l (120071 ). Other proposed features in the density distribution include a rotating 



one-armed spiral density wave based on the variation of the ratio of violet to red int ensities 
in asymmetric double-peaked l ine pr ofiles (see, for example, IWisniewski et al.l 120071 ). For 6 
Sco specifically, ICarciofi et al.l (120061 ) suggest that the disk geometry must change over time 
or that material from the stellar photosphere must be ejected at higher stellar latitudes to 
explain variations in brightness of the system. Models which include these types of density 
enhancements or depletions will necessarily include additional free parameters which will 
require further monitoring of S Sco with contemporaneous observations. In the future, we 
plan to obtain other observables of 6 Sco in order to place further constraints on the radial 
density distribution in the circumstellar disk of the star and to trace the evolution of the 
disk. 
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research. All data and documentation are considered part of the public domain and are 
furnished "as is." These data were collected with the intent of using them for meteorological 
research. The United States government, its instrumentalities, officers, employees, and agents 
make no warranties, express or implied, as to the accuracy, quality, usefulness of the data, 
products, and/or documentation for any purpose. They cannot be held responsible for any 
circumstances resulting from their use, unavailability, or possible inaccuracy. Earth System 
Research Laboratory (ESRL) will make data and/or display products available in World 
Meteorological Organization (WMO), ESRL, or other standard formats which may change 
at times. ESRL reserves the right to suspend or discontinue this service, or portions thereof, 
at any time. 
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Table 1. Standard Stars 



Star 


K 


J-H 


H-K 


Exp. (sec) 


Co-adds 


Gl 406 


6.071 


0.619 


0.364 


0.20 


100 


HD 106965 


7.312 


0.048 


0.015 


0.50 


50 


HD 136754 


7.135 


0.016 


0.005 


0.30 


60 


SA 108-475 


7.986 


0.682 


0.136 


1.00 


20 


HD 161961 


7.304 


0.063 


-0.011 


1.00 


20 


Gl 748 


6.302 


0.524 


0.243 


0.20 


100 



Table 2. 2006 May 12 Photometric Solution 



Extinction Color Term 



Filter ^ points Ext. Ext. Factor Color term Color 

(1) (2) (3) (4) (5) (6) (7) 



J 31 -0.12 ± 0.03 -5.06 ± 0.05 5.26 -0.003 ± 0.023 5.057 ± 0.010 

H 25 -0.16 ± 0.02 -4.31 ± 0.02 2.76 0.004 ± 0.014 4.307 ± 0.006 

K 26 -0.11 ± 0.02 -5.30 ± 0.03 3.11 -0.018 ± 0.017 5.310 ± 0.007 



Note. — Col. 1: Filter. Col. 2: Number of data points used. Col. 3: Extinction 
coefficient (mag airmass^^). Col. 4: Extinction zeropoint and error (mag). Col. 5: Factor 
used to increase errors in Fig. [31 Col. 6: Color term and error (mag color^^). Col. 7: Color 
term zeropoint and error (mag). 
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Table 3. Delta Sco Photometric Solution 



Filter or Color 


Result and Uncertainty 


J 


2.21 ± 0.03 


H 


2.16 ± 0.03 


K 


2.07 ± 0.03 


J-H 


0.05 ± 0.042 


H-K 


0.09 ± 0.042 
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Fig. 1. — Greyscale defoucsed K band image of 6 Sco obtained on 2006 May 12. Although 
the maximum DN is greater than 130,000 the greyscale range is limited to 5,000 to 25,000DN 
to emphasize the donut-shaped image. The full width of the outer ring measured on several 
images is 38 pixels or 14" . The central peak is known as Poisson's spot and arises from the 
diffraction pattern from the secondary. Pixel values are relative to the array center. 
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Fig. 2. — Curves of growth in the H filter for the standard star HD 161961. Left panel is 
before the 5% correction described in §2 was applied to the lower half of the array. Right 
panel shows the same but after the 5% correction. Error bars are comparable to the size of 
the points and are determined from the photon statistics. Different symbols correspond to 
one of the 4 dithered images. 
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Fig. 3. — Plots of fits for extinction (top=J, middle=H, bottom=K) for the standard stars 
observed on 2006 May 12. The mag nitude differerices on the y-axes are the instrumental 
magnitude minus the magnitude from lGuetter et al.l (120031 ) in the indicated filter. These fits 
exclude points which lay within the cosmetically poor regions discussed in §2. Error bars 
have been increased from the photometric error value as discussed in §2 and listed in Table [2l 
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Fig. 4. — Plots of magnitude differences of Fig. [3] corrected for extinction using Columns 3 
and 4 of Table [2] as a function of color of the standard stars. Lines are plotted using color 
term results from Table [2l 
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Fig. 5. — The integrated precipitable water (IPW) content measured at the Camp Navajo 
Army Depot ground based GPS receiver on 2006 May 12. Please note the disclaimer at the 
end of this paper. 
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UT (hours) 

Fig. 6. — The de-extincted J magnitudes of Fig. H] as a function of UT time. Error bars are 
excluded for clarity. The solid line shows the average of the de-extincted magnitudes versus 
time to guide the eye for any trend correlated with Fig. [5l 
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Fig. 7. — The airmass of the observed de-extincted J magnitudes of Fig. [6] as a function of 
UT time. 
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Fig. 8. — Using the photometric solution for extinction in Table [2] we determined the indi- 
vidual S Sco J, H, and K magnitudes. These are shown above as a function of airmass. The 
averages of these points appear in our final solution for 6 Sco in Table [31 
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Fig. 9. — Comparison of the observed Ha lines taken on 19 April 2006 (dotted line with 
circles) and 22 May 2006 (dotted line with squares) and the theoretical lines produced by the 
model with disk parameters n = 4.0 and po =5.0 x 10^^° g cm~^ and with a inclination of 
i = 30° (solid line) and i = 20° (dashed solid line). We also show a model with an evacuated 
annulus between 1.5 and 4.0 stellar radii which had a density of 0.01 compared to our single 
power-law in this region (dot-dashed line). 
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Fig. 10. — Color-color diagram for the constructed models. The single point represents the 
combined photometric observations described in section 2. Each curve represents a different 
power-law density exponent with disk-base densities decreasing from left to right. The star 
shows the best-fit disk parameters obtained by constraining the solutions with spectroscopy. 
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Fig. 11. — Temperature distribution of the disk with model parameters n = 4.0 and po =5.0 
X 10^^° g cm^^. We obtain a minumum temperature of 6700 K in the cool, dense region 
along the equatorial plane of the star and an average overall disk temperature of 12000 K. 
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Fig. 12. — JHK colours from lCarciofi et al.l (120061 ) plotted against the observation described 
in section 2 (point with error bars) and one of the model curves from Figure 9 for a fixed 
n = 4.0 and a disk-base density varying from po = 1.0 x 10^^^ g cm~^ to po = 1.0 x 10^^ 
g cm~^ (dashed line). The JHK colours are grouped together based on time of observation 
(ranging from 2000-2005) and location. The earliest observations are the diamonds, followed 
by the circles, the stars, the triangles and the squares. 
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Fig. 13. — Spectral energy distribution for the secondary star, the isolated primary star and 
the primary star surrounded by a circumstellar disk with parameters n = 4.0 and po =5.0 x 
10-^° g cm-^ 



